Abstract. We present X-ray observations of PSR 0540-69.3 
with the XMM-Newton observatory. The spectra obtained with 
the Reflection Grating Spectrometer reveal, for the first time, 
emission from ionized species of O, Ne and Fe originating from 
the SNR shell. Analysis of the emission line spectrum allows 
us to derive estimates of the temperature, ionization timescale, 
abundances, location, and velocity of the emitting gas. 
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1. Introduction 



2. Observations & Data Reduction 



SNR 0540-69.3 is one of only a few Crab-like supernova rem- 
nants (SNRs) with a shell. SNR0540-69.3 harbors a pulsar, 
PSR B0540-69.3, disc overe d in the soft X-ray band by Seward, 
Harnden, & Helfand ( 1984 ). Like the Crab, it shows plerion 
emission from a synchrotron nebula powered by the embed- 
ded young pulsar (Chanan, Helfand & Reynolds 1984). From 



the period (50 ms) and spin down rate (4.79 x 10 ss ) a 
a characteristic age of 1660 yr and rotational energy loss rate 
of 1.5xl0 38 ergs s _1 have been derived (Seward, Harnden, & 



Helfand 1984) 



Early optical observations (Mathewson et al. 



1980J) classi- 
fied SNR 0540-69.3 as a young oxygen-rich SNR with an 8" 



diameter shell, bright in [Olll]. Mathewson et al. (198C) also 
found a diffuse patch of [O III] ^30" to the west of the ring 
which indicated that the remnant was larger than 8". A ROSAT 
High Resolution Imager observation by Seward & Harnden 
(1994) also revealed emission from well outside the central re- 
gion, which they interpreted as from a patchy outer shell with a 
diameter of 55". This shell contributes ^20% to the measured 
flux in the R OSAT 0. 1-2 keV band. More recently Gotthelf & 
Wang ( 200C ) presented a high resolution Chandra HRC obser- 
vation of PSR B0540-69.3 which clearly shows emission from 
the outer shell. 

In this letter, we report the results of observations of SNR 
0540-69 acquired with the Reflection Grating Spectrom eter o n 
board the XMM-Newton Observatory (den Herder et al. 2001 ). 
These observations allow us for the first time to detect and iden- 
tify emission lines from the remnant shell. Our spectral analy- 
sis allows us to derive values for the temperature, ionization 
timescale, abundances, location, and velocity of the emitting 
gas. 
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* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and the USA (NASA). 



PSR B0540 was observed repeatedly as part of the calibration 
of the Reflection Grating Spectrometer (RGS) during the Cal- 
ibration/Performance Verification (Cal/PV) phase. For the pur- 
pose of this letter, we have analyzed the 126 ks data obtained 
on 2000 May 26, which is the longest observation with the low- 
est background. The observation was performed with the tele- 
scope rolled such that the RGS dispersion axis was aligned at 
21.29°, clockwise on the sky, from celestial north. The data 
were initially processed with the XMM-Newton Science Anal- 
ysis Software (SAS). The spectra were extracted by applying 
spatial filters to the CCD image while a CCD pulse-height filter 
was applied to select the m = — 1 spectral order. The remnant 
is located on the edge of a broad region of diffuse emission 
(Wang & Helfand 1991 ) which complicates background deter- 



mination. However, the image is small enough on the detector 
that spatial regions offset from the source position can be used 
for background subtraction. The background spectra were ex- 
tracted by applying the same pulse height filters as used for the 
source, to spatial regions at either edge of the camera. 

A response matrix appropriate to the spatial extent of the 
source was generated as follows. A spatial mask corresponding 
to the RGS aperture was imposed on the Chandra HRC im- 
age, and the intensity distribution was integrated over the RGS 
cross-dispersion direction. The resulting profile was convolved 
with the RGS point source response matrix, generated with the 
SAS task RGSRMFGEN. The spectral analys is was performed 
using the SRON SPEX (Kaastra et al. |1996|) p ackage, which 
contains the MEKAL code (Mewe et al. 
thermal emission. 



3. Analysis and Results 

3.1. Location of emitting region 

We initially extracted spectra with spatial filters 0.25', 0.5', 1' 
and 1.5' wide, each centered about the peak emission. The most 
interesting feature in the spectra is the presence of an O VIII 
Lya emission line that becomes more prominent as the width 
of the spatial extraction region is increased. This is an indica- 
tion that at least some of the O emission is from the SNR shell, 



1995) for modeling 
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Fig. 1. Cross-dispersion profiles of SNR 0540-69 obtained by 
integrating along the dispersion axis for a: the entire spectrum 
and b: O VIII Lya only. 



and not the central region. To further constrain the spatial dis- 
tribution of the emitting gas we created a cross-dispersion im- 
age profile of the O VIII Lya line by integrating over the line 
in the dispersion direction. The cross-dispersion profile, plot- 
ted in Fig. [lj, shows two peaks. Referring to the Chandra HRC 
image (Gotthelf & Wang 2000} reproduced here as Fig. 2), we 
conclude that we see O emission both from the central region 
of the remnant, centered on the pulsar, (the peak at offset zero 
in Fig. lb), as well as from the partial shell to the west (the 
second peak sa 24 arcsec towards more negative offset in Fig. 
lb). 

The angular extent of this partial shell, A(j> w 30", corre- 
sponds to an effective wavelength width along the dispersion 
direction of * 0.06 A, comparable to the wavelength resolu- 
tion of the spectrometer. This implies that we do not have the 
sensitivity to resolve variations in brightness along the shell. 



3.2. Spectral analysis 

We divided the remnant into three regions, Northeast (NE), 
Central (C) and Southwest (SW), for the purpose of spec- 
tral analysis. For the NE and SW regions we use an extrac- 
tion width of 0.5'centered at ±0.375' about the central com- 
pact region, while for the brighter C region we use a width of 
0.25'. The spectra for these regions are displayed in Figs. 3 
and 4. Each spectrum shows a strong power law continuum 
with superimposed line emission, most notably Ovill Lya. 
Emission line features are weakest in the spectrum extracted 
from the NE region, where we only detect O VIII Lya and pos- 
sibly weak Ovil emission. The line emission is strongest in 
the SW spectrum where we easily identify emission lines of 
O VIII (Lya, Ly/3), O VII (n = 2-1 triplet, He/3), Fe XVII 
(15.01,16.78,17.05,17.10 A), Nex (Lya), and Neix (n = 
2 — 1 triplet). Due to galactic absorption, no features can be 
seen longward of ~ 23 A. Note that even though the shell- 
like structure and the plerion are only partially resolved by the 
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Fig. 2. X-ray image of SNR 0540-69.3 obtained by the Chan- 
dra HRC. The image contains the 50 ms pulsar PSR B0540- 
69.3, the plerion and the outer shell. The SW, C and NE regions 
used for spectral analysis are indicated. The disperion direc- 
tion of the RGS is along the long dimension of the extraction 
regions. 



telescopes of XMM-Newton, the high spectral resolution of the 
RGS allows for high sensitivity to the discrete emission from 
the remnant, even though the spectrum is dominated by the 
powerlaw continuum from the plerion. 

In the SW spectrum, O VIII Lya has a significant offset of 
AA = —0.17 A with respect to the laboratory wavelength. This 
offset must be largely due to a Doppler shift, and not to a posi- 
tional offset between the nominal pointing direction (the pulsar 
position) and the emission centroid of the shell, because the 
equivalent wavelength extent of the shell is much smaller than 
the observed shift. In the Chandra image, we find that the cen- 
troid of the S W emission, projected onto the RGS dispersion di- 
rection, is offset from the centroid of the plerionic emission, by 
about w —9" along the dispersion direction. This corresponds 
to a wavelength offset of —0.019 A. The net Doppler shift of 
the SW emission is thus —0.15 A, which corresponds to a ra- 
dial velocity of v ks —2370 km s _1 . A similar radial velocity 
is also observed in the other strong emission lines in the spec- 
trum. Note that the radial velocity of the LMC, wlmc = +278 
km s _1 , is only a relatively small correction to this measure- 
ment. The emission lines seen in the C and NE spectra are too 
faint to permit a meaningful radial velocity determination. 

We synthesized a spectral model comprising three compo- 
nents: a power law to represent the synchrotron emission from 
the plerion, a non-equilibrium ionization (NEI) model for the 
thermal line emission, and foreground absorption. The free pa- 
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Fig. 3. RGS 1 & RGS 2 (upper and lower plots respectively) 
spectrum of the SW region. The solid line is a model fit to 
the data. The gap between 20 to 24 A in the RGS 2 spectra 
is due to the failed CCD 4 on that instrument. The spectra ex- 
hibit a strong power law continuum due to synchrotron radi- 
ation from the plerion and thermal line emission from highly 
ionized atoms of O, Ne and Fe 



rameters of the model are the photon index of the power law 
and its normalization, the electron temperature T c , emission 
measure, abundances, ionization age n c t of the shocked gas, 
and the column density TVh of foreground absorbing gas. Here, 
n e is the electron density of the shell, and t is the time since the 
hot gas was shocked. We adopt a distance of 5 1 kpc throughout. 

We concentrated our quantitative analysis on the S W region 
since the spectrum from this region contains the strongest line 
emission features. Fixing the O abundance to a third solar (note 
that we have no independent way of constraining the absolute 
abundances, since the continuum emission is dominated by the 
plerionic emission), the abundances of Ne and Fe were allowed 
to vary with respect to O. An initial fit was made to the SW 
spectrum, using O VIII Lya and O Vll-triplet only. The param- 



eters obtained were subsequently used as starting parameters 
for fits to the entire spectrum. For the Ne/O and Fe/O abun- 
dance ratios, we obtain Ne/O ~ 3.0 and Fe/O ~ 2.1 relative 
to their solar values. The best fit model parameters are listed 
in Table 1 . Fits were also made to the NE and C spectra. Since 
the line emission in these regions is extremely weak we fixed 
the NEI component parameters (i.e. electron temperature and 
ionization parameter) to the best fit parameters obtained from 
the SW region, while allowing the emission measure and pho- 
ton indices of the power law component to vary. The best fit 
parameters are given in Table 1, and the spectra, together with 
the best fit models are shown in Figures 3, 4 and 5. 

As can be seen from the Figures and the \ 2 values in Table 
1, the fits are not perfect. This is probably due to a combination 
of remaining small calibration uncertainties in the RGS effec- 
tive area (most of the \ 2 is m the relatively poor fit to the con- 
tinuum, which is due to the plerion, at the shorter wavelengths, 
where the RGS effective area is least well calibrated), and the 
simplicity of our spectral model, which does not allow for tem- 
perature and ionization inhomogeneities. In particular, there is 
a feature at w 14 A, which is not fit by our isothermal model. 
Candidate emission lines are 2p — 3d in Fe XVIII at 14.208 A, 
or possibly 2s - 3p in Fe XXI at 14.008 A, although the latter 
should be accompanied by strong emission at 12.29 A in equi- 
librium. In any case, detectable amounts of higher charge states 
of Fe L would definitely indicate the presence of hotter plasma. 



Table 1. Fitting results for the SW, C and NE regions of SNR 
0540-69. 



parameter 


SW 


C 


NE 


NEI: 

n c n H V (10 58 cm- 3 ) 

kT e (keV) 

n c t (10 9 cm" 3 s) 


92±6.5 

0.58±0.18 

25.22±5 


50±20 


45±18 


POW: 

photon index 


1.82±0.1 


1.2±0.04 


1.75±0.1 


ABS: 

N H (10 21 cm 2 ) 


3.88±1 






X 2 /d.o.f. 


2.04 


2.8 


2.9 



4. Discussion 

First of all, a qualitative point. We detect faint, but significant 
O VIII emission from the direction of the plerion. At our spa- 
tial resolution, we cannot exclude the possibility that a signifi- 
cant component of this line emission is in fact from foreground 
remnant gas seen in projection and contamination by the SW 
region. If confirmed at higher spatial resolution, this detection 



would be interesting, because Mathewson et al. ( 1 980 ) detected 
intense optical [O III] emission from this same inner region, 
with a large velocity dispersion of v « 3000 km s _1 . There 
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may consequently be a wide range of ionization present in the 
plerion. Unfortunately, we do not have the sensitivity to try and 
localize the highly ionized oxygen, and to constrain its kine- 
matics for comparison with the [O III] image. 

Next, we check for consistency of our measured parame- 
ter values. The interpretation of the radial velocity measured 
in O VIII is somewhat complicated. We see what appears to be 
a limb-brightened shell in the Chandra HRC image, but if the 
emission is really arranged in such a shell, one would not ex- 
pect to see any radial velocity shift of the material at all, its 
velocity vector being entirely perpendicular to the line of sight. 
We know instead that the line emitting material has a velocity 
component along the line of sight of w —2370 km s _1 , which, 
combined with an unknown perpendicular velocity component, 
sets a lower limit on the true space velocity of the O VIII emit- 
ting gas. To relate this space velocity to the velocity of the ex- 
panding blast wave requires making a further assumption about 
the structure and evolution of the expanding remnant. However, 
the measured radial velocity still gives a lower limit on the true 
expansion velocity, independent of these assumptions. 



With a measured radius of s=a 30", the shell has a linear 
radius of R s = 7.4 pc. Combined with the lower limit on the 
expansion velocity, we estimate an upper limit to the age of the 
shell of t < 2800 yr, consistent with the pulsar spindown age 
of 1660 yr. Conversely, if we assume that the age of the shell is 
identical to the pulsar spindown age, we would conclude that 
the true expansion velocity of the blast wave is w 4600 km 
s , and the radial velocity we observe in the X-ray lines is 
either due to material well inside the interior of the remnant 
(where the flow velocities are smaller), or the line emitting gas 
must have a substantial perpendicular velocity component if it 
is located in the immediate post-shock region. 

The discrete emission line spectrum implies an ionization 
age of approximately n c t ~ 2.5 x 10 10 cm~ 3 s, a low value, 
which probably implies that the plasma has not yet reached ion- 
ization equilibrium. Combined with the pulsar spindown age, 
one infers a density of the medium of n c ~ 0.4 cm~ 3 . But 
the gas may have been shocked (much) more recently than the 
characteristic age of the blast wave, depending on it position 
within the remnant interior, and so we should regard this den- 
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sity estimate as a lower limit. Also, from the emission measure Seward F.D, & Harnden Jr F.R., 1994, ApJ, 421, 581 
of the line emitting gas, we derive a density of approximately Wang Q.D., & Helfand D.J., 1991, ApJ, 379, 327 
n c ~ 10 cm~ 3 (assuming a half shell of relative thickness 
AR ~ i? s /10 for the volume estimate). Given the relative in- 
sensitivity of the density estimate to the volume and emission 
measure, this higher density estimate may be closer to the char- 
acteristic values at the shell than the estimate derived from the 
ionization age, and this would favor the conclusion that the gas 
has been shocked more recently. In fact, the measured electron 
temperature, kT c ~ 0.6 keV, is incompatible with the lower 
limit on the blast wave velocity v s > 2400 km s _1 (unless the 
electron- and ion temperatures have not equilibrated), and most 
likely, we are seeing material that has recently passed through 
a (much slower) reverse shock. If gas with a much higher elec- 
tron temperature is associated with the outer blast wave, it will 
be difficult to detect it in the XMM-Newton data. Because the 
remnant is only partially resolved, both the RGS and EPIC data 
are dominated by strong emission from the plerion, which re- 
duces the sensitivity to detection of hot thermal gas. 

5. Conclusions 

We have detected resolved X-ray line emission from the su- 
pernova remnant shell around PSR B0540-69.3, confirming 
the nature of the shell-like structure seen in previous exper- 
iments as a thermal SNR. Unfortunately, given the relatively 
low signal-to-noise of the observation, the fact that the object 
is only partially spatially resolved, and the uncertainties in the 
interpretation of the estimates of the various observable param- 
eters, we can not uniquely decide what the evolutionary status 
of the remnant is, or what the relation of the line emitting gas 
is to the overall structure of the remnant. 
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